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riu'  .id.iptivc  Spi'ilral  l.ino  I.  ii  h .1  lu  o r (Sl.l  ) is  a ilovii.f 
f 0 r a II  t oma  t i c a 1 1 >■  p rc-w  h 1 t cn  i ng  noise  uifh  unknown  speetriim 
a lul  extracting  (enhancing)  spectral  comjioncnts  to  which  the 
noise  is  additive.  The  steady  state  nerformance  of  the  Sl.l. 
presented  for  the  general  case  of  multiple  s [) e c t r a 1 lines. 
I'he  theoretical  steady  state  S.NK  performance  of  the  Sl.h  is 
sliown  to  be  that  of  a sjiectrum  analyzer  with  resolution 
determined  b >•  the  1 e n g t )i  o f the  a d a p t i v e 1 >'  weighted  t a p p e d 
delay  line  which  is  used  to  implement  the  tlevice.  _ 
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I ()  I-  s i n a 1 tl  i- 1 oc  t i o n prnlilciii'.  of  int<.'rcst  in 

onar,  i dotoction  processor  is  lequired  to  detect  sjiec- 
tral  lines  in  the  presence  of  an  add  i t i v i'  broadbaiul  noise 
b a c k g r oiin  d . The  optiinum  { I.MS  ) detector  in  this  instance 
requires  exact  knowledge  of  both  the  signal  a lul  noise 
spectra.  However,  in  general  this  type  of  a priori  know- 
ledge of  signal  and  noise  spectra  is  not  available.  If 
the  signal  is  b a n d 1 i m i t e d ( i . e . spectral)  a n tl  the  noise 
IS  broadband,  then  the  signal  c o m ]i  o n e n t is  c o r r e 1 1 a t e d 
for  time  lags  which  are  large  compared  to  that  for  the 
noise.  This  difference  in  correlation  time  can  be  exploit- 
ed to  build  a self-adjusting  filter  which  will  automatically 
extract  the  spectral  c o m |i  o n e n t from  the  b r o a d b a n d noise 
i n d e ji  e n d e n t of  the  frequency  of  the  spectral  c o m j)  o n e n t . 

Such  an  adaptive  filter  has  been  called  ;i  "Spectral  |,ine 
I^nhancer  (Sl.li)  [I)-"  This  document  presents  a steady 
state  analysis  of  the  S 1. 1;  , derives  the  adaptive  filter 
control  algorithms  and  convergence  characteristics  and 
considers  the  limiting  case  of  a white  (iaussian  noise 
nVG.N  ) background. 


The  intent  of  this  document  is  to  show  that  the  Sl.T. 
is  an  a d a [1 1 i V e noise  p r e - w h i t e n e r and  signal  notch  filter 
which  automatically  follows  any  spectral  components  in  the 
]i  r e s e n c e of  broadband  noise.  It  is  observed  that,  in 
fact,  the  local  ,S N' R at  the  output  of  the  notch  filter,  i . e . 
within  the  band  defined  by  the  resolution  of  the  ada]itive 
filter,  is  not  increased  but  rather  remains  const, int. 
finally,  in  general,  it  is  intended  to  give  s o nu-  additional 
insight  into  the  .S  I.  )•;  operation  to  assist  both  in  exploiting 
its  potential  and  understanding  its  limitations. 


II.  STHADY  STATh  ANALYSIS 


A func  t i ona 1 
( S 1. 1;  ) is  given  in 


diagram  of  the 
I''  i g 11  r e I . T h e 


s p e c 1 1 a 1 line  e n h a n c e r 
signal  ['Ills  noise  input 


X ( t ) = s ( t 1 + n ( t ) 


( 1 ' 


I s p roc  ess  ed  i n 
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to  il  o - c o r r c 1 a t o tho  broadband  additivo  iioIm'  component 
n^(t)  from  n (t  A),  and  an  adaptively  wei^;htcc.l  l.ipped 
delay  line  to  |)rovide  a minimum  Mean  .Sqii.ired  Irror-  ( '!  S I i 
estimate,  s ( t ) of  the  s i n a 1 s f t ) . The  n j y,  n ;i  1 ■ a n c e 1 1 a 

t i on  channel  subtracts  the  MSH  estimator  s(tj  from  the 
s i n n a 1 - p 1 u s - n o i s e input  to  provide  an  error  s i y n . i 1 ( t J . 

fhe  power  in  the  error  sij’nal  is  then  minimized  with  res- 
pect to  the  weiuhts  in  the  tapped  delay  line.  Since  the 
noise  in  the  s i j;  n a 1 c ;i  n c e 1 1 a t i o n c h ;i  n n e 1 , n f t J , is  uncor- 
re  lilted  with  the  noise  in  the  signal  estimation  channel 
n 1 ( t ) = n ^ ( t - 6 ) , the  process  of  minimizing  the  error  out 

put  power  works  primarily  on  rediicitiK  the  sinn.al  power  in 
the  error  output.  The  MSI;  estimator,  s(tj,  is  the  desired 
Sl.l:  output.  this  estimator  represents  a " c 1 ea  n ed  - ii  p " , i.e. 
low  additive  noise,  version  of  the  input  siyiial  i.ompoiient 
s ( t ) which  will  be  cancelled  b >•  subtracting  s f t ) . 


Let  the  signal  component  of  the  input  consist  of  '■I 
spectral  terms  and  be  represented  b \' 

s ( t ) = }"k.  (t  )e-"^i^ 

i = 1 ‘ 


C) 


where  Sj(t1  is  the  comple.\  envelope  of  .a  spi'etr;il  sign  a 
at  frequency  oj.  . If  the  s u[ie  r s c r i [i  t "1"  is  the  m.itrix 
transpose  operator,  then  the  signal  plus  noise  vi'ctor  X 
is  defined  where 


=[.x(t-A)x(t-6-A) 


X ( t - ■ 


■)] 


(.M 


= sT  a n'I' 


(t) 


such  that 


= fs(t-A)s(t-6-A)  •••  s(t-A-[l,-lJAjJ 


(5) 


= [n^^  ( t - 6)  n ( t - 6 -Aj  •••  n ( t - 6 - ( L - 1 ] A )].  (6) 


A signal  delay  matrix 

^ '**  UmI 


(7') 
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in  coniunction  with  a si^;nal  complox  cnvcloin-  vector  liefineil  a 


I t 

f 1 ) ■■  ' ? • • • ■ 1 ' I 


I M'j 


|M  r III  I I r ew  I I t I II  i;  tin  I ?■  ,i  ii  |i  n •,  e n I ( 1 ) a 

X III.  ‘ N. 


0 I 


0) 


Notice  that 


s ( t ) = 1 


(11) 


where  1 ^ = 1 1 1 " ' * 1 1 is  an  M - e 1 e m e n t vector.  1 1 e f i n i n g the 

mat  r i x'comp  1 ex  conjugate  transpose  operator  "II",  a weight 
vector  W IS  given  by 

w"  - 

[ w j * w 2 * • • • w [ * 1 

(12) 

where  " * " denotes 
is  now  expressed 

the  complex  conjugate, 
by 

rhe  error  output 

f(t)  = 

x(t  ) - Sft) 

( 1 .3  1 

= 

X f t ) - w'*l)  h - w"n 

Ml) 

with  expected,  i . 
signal  and  noise 

e . steady  state  power, 
given  by 

f o r u II  c o r r e 1 a t ed 

= f:  M ) I ‘ ^ - 2Re^w"l)t(H  - 


( 1 la  ) 


In  (II),  the  fol  lowing  definitions  have  been  made 


E<  • 1 

= s t ;i  t 

i s t i c a 1 

c X p e c t a t i o n 

o]ie  r a t o r 

(15  1 

EUt"} 

= R 

s s 

( s i g n ;i  1 

c o r r c 1 a t ion 

m ;i  t r i X ) 

( 1 bl 

t { NN** } 

= Rvv 
NN 

(noise 

col  elation 

m ;i  t 1'  i X ) . 

( 1 ') 

TM 
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1 1 1 DO 

The  minimum  va  1 ue  o 1' 
obtained  by  solving 

L M ' ( M ! " 1 

1 with  re''  pe  e 1 

! t o 

W 1 

{ 1 '.  ( t ) r ) 

9W 

= -2DR 

s ^ 

= 0 

1 + 2|Dli  d" 
; s s 

* R 

1 W 

NN  - 

I 


f 1 « ) 


for  W to  yield 


Nisr.  = '""ss"  ^ *'nn1  '''Ssi- 


f I 9 ) 
(20) 


Using  the  identity 

(m**q"*m  + p'* 


H II  . . 1 

P - PM  [MPM  *•  Q|  MP 


(2  1 } 


(2G)  becomes 
. 1 


^ISE  " '^NN 


1)(  I 


H - 1 ■ 1 

[UR  'd  + R 
' NN  ss 


- 1 


DR  D)R 
N N ^ ^ ■ 


(22) 


To  obtain  an  expression  for  the  M S E , E ( | r ( t ) | “ } , and  signal 
estimator  channel  output  power, (22)  is  used  in  (14)  aiui 


s 1 t ) 1 } 

= 6f 

Iw"  (D1-;  + 
-MSP.  - 

N)  1 ^ t 

( 2 .-S  ) 

res 

pec  t 

i ve 

1 y to  y 

i e 1 d 

. 1 

E'  1 

e (t) 

!^) 

MSP 

S + N 

- 1 

[DR^sD"  ^ R^^^^l 

(24  ) 

and 

'dNl'^”'dsi- 

Ff  1 

y (t  ) 

1^} 

MSIi 

T 

1 «ss 

D[I)R 

[)**  + 
s s 

(2S  ) 

1 n 

(24) 

the  s i g n a 

1 and 

no  i 

sc  p 0 w 

er  have  been  defined  as 

S = 

Eds 

(t)  1 

d 

( 2b) 

and 

N = 

E(  In 

o('> 

7 

(27) 
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t<|uation.s  (22)  and  (25)  arc  the  jjciici  il  expressions  for 
the  optimum  ('ISli)  Sl.f  filter  weiRhts  and  expecteii  Sl.f  output 
power  r e s p e c t i V e 1 V for  M spectral  components  in  ,a  correlated 
noise  background.  The  following  special  cases  will  be  con- 
sidered; (a)  small  signal  to  noise  ratio  and  (b)  a single 
spectral  component.  The  small  signal  to  noise  ratio  case 
will  indicate  that  the  general  structure  of  the  Sl.f,  signal 
estimation  channel  is  the  classical  "(noise)  pre-whiteri  and 
(signal)  match"  filter  operator.  The  single  spectral  com- 
ponent case  will  illustrate  that  the  .Sl,li  signal  matched 
filter  function  is  implemented  by  a simple  Pi  sc  ret  e I'ourier 
Transform  (DFT)  filter  with  effective  5dB  filter  width,  i.e. 
resolution,  equal  to  (1/1. A)  Mz  and  center  frequency  e(|ual  to 
the  frequency  of  the  spectral  comiionent.  Recall  1.  is  the 
number  of  filter  taps  and  A is  the  temporal  sjiacing  between 
taps.  Small  signal  to  noise  ratio;  As  the  signal  power 
terms  on  the  diagonal  of  R^^  becomes  small  relative  to  the 
noise  power  terms  on  the  diagonal  of  R«g^^■.  the  MSI.  weight 
vector  I 2 0 ) and  S 1,  t expected  output  p o w e r £ M y ( t ) ' “ .‘  (25) 

approach 


and 


—MSI-.  N N s s — 


(2R  1 


(29) 


First,  the  structure  of  the  S 1. L signal  estimation  c h ;i n n e 1 
as  dictated  by  (28)  is  examined.  Figure  2 gives  a functional 
block  diagram  of  the  SLh  processor  and  indicates  the  impor- 
tant filter  o ]i  e r a t i o n s in  the  signal  estimation  channel. 
Observe  that  the  matrix  filter  operation,  R ' . , on  the  delayed 
data  vector  X is  clearly  a noise  ['  r e w h i t e n i ti  g function.  T h e 
scalar  filter  operation  is  equivalent  to  a signal  matched 
filter.  This  can  be  seen  bv  examining  the  SI. I output  in 
more  detail.  The  output,  s^ ) , is  given  by 


T M H - 1 

s(t)  = i'r  "d  R^‘x 
— s s NN  — 

= I P I’nX 
^ Inm— n-Tv 
n , m = 1 


(30) 
(3  1) 
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whore 


nm 


= 


(t  ) s_ 


(t ) } 


j ((‘'n  -<‘'m  ) 


is  the  correlation  between  the  and 
ponent,  complex  envelopes  and  = R*' 
as  the  noise  prewhitened  data  vector. 


f ) 

1 1)  , 
n spectral  c o m - 

^ is  referred  to 

The  term 


a 

n 


= n«x 

— n — w 


Lo  w 


(t-  [£-l  ]A)e 


£=  1 


1 6 
‘^n^e^  n 


(.^3) 

(34) 

(35) 


(XI  j-» 

since  X , , is  recognized  as  the  Discrete  ^ourier  Transform 
(DF'T)  oY  the  noise  prewhitened  data  waveform  x^^{t).  The  OFT 
is  evaluated  at  frequency  As  indicated  by  (31),  the 

output  of  the  n^^  DFT  filter  channel  is  weighted  b>-  the 
factor  M and  summed  together  with  the  outputs  of  the 

^nm 
M=  1 

(M-1)  other  correspondingly  weighted  DFT  filter  outputs. 
Notice  that  for  correlated  tonals,  i.e.  0 for  n^m, 

coupling  between  DFT  filter  outputs  occurs.  This  umlcsirc- 
able  effect  could  lead  to  time  varying  (intermodulation) 
effects  characterized  by  periodic  signal  suppression  due  to 
destructive  combining  of  out  of  phase  spectral  components. 
This  type  of  behavior  is  a distinct  possibility  in  an  adap* 
tive  realization  of  the  .SLF  wherein  the  requirement  for  short 
adaptation  time  for  realization  of  see  Section  IV) 

allows  for  residual  correlation  between  complex  envelope 
functions  which,  in  fact,  become  more  uncorrelated  only  as 
the  correlation  averaging  time  increases.  The  ideal  situa- 
tion, of  course,  is  where  D =0  for  m^n  such  that 


s (t) 


M 

= I X 

^ , m m — m — w 
1 


(36) 


and  each  DFT  channel  output  is  weiglitcd  proportionately  to 
the  power  in  the  spectral  line. 
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The  SLIi  output  jiower  E{|y(t)l  ’ the  small  siKnal 

to  noise  ratio  case  with  uncorrelatei.1  sj)cctral  component 
envelopes  and  uncorrelated  noise  can  be  obtained  from  (29). 
For  this  case  R is  a diajjonal  matrix  with  m**^  element 
p , the  OFT  matrix  is  such  that  = <5  fdirac  delta), 

il'e.  orthogonal  carriers  and  ^NT  whc?P  I is  an  I.  by  I, 

identity  matrix.  Accordingly,  (29)  reduces  to 


carriers  and  Rj^  = ! 

Accordingly,  (29) 
M 


m element 
fdirac  delta) 
I is  an  I.  bv 


F.{|y(t)|'}|M,,  = 

It  is  important  notice  that  according  to  (.^7)  the  only  way 
to  increase  the  SLF.  output  signal  to  noise  ratio  is  to 
increase  the  resolution  of  the  adaptive  ni-T  filters,  i.e. 
by  increasing  the  length  of  the  adaptive  tapped  delay  line 
L in  F igure  1 . 

An  exact  expression  for  SLF;  output  power  for  a single 
spectral  line  in  additive  uncorrelated  noise  is  obtained 
from  (25),  namely  (p , , = S) 


namely 


E{|y(t)r} 


Therefore,  even  for  small  S/N  if  the  tapped  delay  line  length 
L is  made  sufficiently  large  the  spectral  line  can  be 
"enhanced"  to  the  extent  allowed  by  the  resolution  of  the  1.- 
point  PFT.  The  MSB  weight  vector  case  of  a 

single  spectral  line  in  uncor  re  1 at  ea'^rio  i se  is  obtained  from 
(20)  as 

^ f'  (-^9) 

-MSfc  N I I q ^ 


By  noting  that  £{|y(t)| 


of  the  form 


£{  |y(t)  I } = - /3^LN, 

the  SLE  output  signal  to  noise  ratio  is  given  by 

SNR  = tf 
out  N 


= L^NRj 


-4 
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Thus,  for  example,  a single  spectral  line  with 

can  be  enhanced  to  = 0 dB  with  an  I,  = 100  tap  delay 

line. 

A qualitative  example  of  precisc'ly  how  the  Sl.T  would 
function  for  a hypothetical  spectral  signal  and  broadband 
noise  is  informative.  f'igures  3 f a 1 , (b)  and  (c)  give  the 

power  spectrum  of  the  SLF.  input,  the  noise  prewhitener  out- 
put and  the  matched  filter  output.  A signal  consisting  of 
two  spectral  lines,  plus  "colored"  broadband  noise  is  con- 
sidered in  F-'igure  .3(a).  The  colored  noise  is  c h a r a c 1 e r i c ed 
by  a sloping  spectrum  which,  as  a result  of  the  p r e w h i t e n i n g 
matrix  filter  operation  is  "magically"  transformed  into  a 
flat  noise  power  spectrum  of  level  N (Figure  3(b)).  The 
I) FT  matched  filters  then  center  at  tiie  line  frequencies  fj 
and  f2  and  exhibit  filter  characteristics  determined  by  i' 
3dB  main  lobe  of  1/LA  (Mz)  and  13.6  dB  down  first  side- 
lobes,  i.e.  a uniformly  weighted  1)1  T.  Notice  that  on  a 
per  Hertz  basis  at  f^  and  f2  the  SNR  remains  at  2dB  and  4dB 
respectively.  However,  the  broadband  SNR  theoretical  !>•  con- 
tinues to  increase  as  I.  increases. 
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1 V . Tilt':  Sl.t;  ADAOTIVI:  lUiAl.lZ  Aj;^!  ON 

All  adajitive  least  mean  s(|iiare  fl.MS)  realization  of 
the  SI, I iitilizinj;  the  method  of  .steepest  rj>radientl  des- 
cent attempts  minimize  the  error  output  power  at  time  t; 
with  respect  to  the  weight  vector  W | t = t j = ^ ( i ) . There 
fore, 


Uit^r'  = 

[xftj)  - w"fi)X(i)J  j^x(tj)*  - x"(i 

1 )W(i  )j 

(4. -5) 

= 

txft.)|“  - 2Re( w" ( i ) X f i ) X f t * 1 } + 

1 w"(i ) X (i ) 

! “ • (44) 

The  gradient 

of  ( 4 .^ ) with  respect  to  W ( i ) is 

3 1 E(ti  ) 1 ^ _ 

U({) 

- 2X(i)x(t.)*  2 X ( i ) x"(  i )W(  i ) 

(4  5) 

= 

- 2X(i)t*ft.) 

(4h) 

which  determines  the  steepest  descent  tapped  delay  line 
weight  vector  control  algorithm 


W(i  + 1 ) = W ( i ) - 2UX  f i)  e*(t  . ) . 


(47) 


In  (47),  the  constant  0 determines  the  rate  of  convergence 
to  a steady  state  estimator  and  the  stabilit)'  (adaptation 
noise)  of  the  final  estimate  of  ^ i.  • T t*  e .x p e c t c d value 
of  (47)  after  K iterations  of  (4^'  is 

E(W(K)}  = M(K) 

“nn’] 

shown  to  converge  to  (20)^ 
i.e.  if  u relative  to  the  m.iximum  eigenvalue  \j,  of 

l>K  .P  ♦ Ryv  satisfies  the  bourn! 
s s N N 

1 4 


= [ 


I ♦ 2u  (im  p + 
s s 


K - If 

2u);  I +211  (PRj.^p'* 

k = o*- 


which  in  the  limit  K c a n b e 


It ' 
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- ‘ n ■ 0.  (SO) 

^ 1 

lor  uncorre  1 at  oil  spectral  onvolo|>os,  ii  n c o r r o 1 a t oil  noise 
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The  material  in  this  section  follows  previous  work  (for 
example  [-])  exactly  and  the  reader  is  referred  to  this 
work  for  additional  results  on  con verj; once  properties  and 
steady  state  m i s ad i u s t men t . 


V . SUMMARY 

The  Spectral  I^inc  l^nhancer  is  a device  for  adaptive 
broadband  noise  re  whitening  and  spectral  line  extraction. 

It  has  been  shown  that  because  the  SI.l:  is  essentially  an 

adaptive  discrete  Fourier  Transform  Filter  (UFT)  bank  (one 
I) FT  channel  per  spectral  line)  its  ultimate  steady  state 
SNR  performance  is  no  better  than  that  determined  by  the 
resolution  of  the  DF-'T.  F’erhaps  the  most  important  aspect 
of  the  SLF  is  that  it  performs  a background  noise  norm. a lie  a - 
tion  function  by  virtue  of  its  prewhitening  capability. 
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From:  Comma  rui  i n s Officer 

To:  nistribiition  Fist 

S 11  b i : Technical  M e m o r a n il ii m ; forwarding  of 

F.ncl:  (1)  N'USC  Technical  Memorandum  No.  Tl)  1 1 1 - 1 0 0 - 7 4 , ”.\  Steady 

State  .\na  lysis  of  the  Adaptive  Spectral  Line  linhancer" 

1.  The  enclosed  document  (1)  is  forwarded  for  your  information 
and  retention. 

2.  This  document  entitled  "A  Steady  State  Analysis  of  the 
Adaptive  Spectral  Fine  F.nhanccr  (SFF)",  was  written  as  a follow- 
up to  a discussion  between  .lohn  Neely  (NAVSMIPS  PMS  .'^02-4)  and 
the  author. 

7).  The  intent  of  the  document  is  to  provide  PMS  ,T02-4  with  a 
more  substantial  theoretical  background  of  the  SFF.  hardware 
which  is  under  development  at  NIIC . 
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